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We make statistics about the number distribution of normal radio pulsars according to the pulsar spin period. And find 
the physical meaning of the traditional statistical method about this problem is very limited. According to the statistical 
method proposed by us, the distribution rule is exactly the exponentially decay with the spin period increasing, especially 
important, which can be interpreted as normal radio pulsar number decays as the radiative elements decay. It discloses that 
though belonging to two different extreme scales of the universe and having different inner active mechanisms, the normal 
radio pulsars and radiative elements accord with the same varying rule. And like the half decay period of the radiative 
elements we can also introduce the half decay spin period of the normal radio pulsars. 
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1 Introduction 

Pulsars — a kind of charming compact objects in the sky, 
are generally known to be born in the supernova explo- 
ration. Exploring them will provide the opportunity for peo- 
ple to understand the fundamental physics and astrophysics 
(Cordes 2004). There are many questions about pulsars de- 
served to be studied. In this paper we would concentrate on 
one of the attractive problems that is how the pulsar number 
distributes according to the spin period. In section 2 we in- 
troduce the usual statistical method at present on the spin 
periods of normal pulsars; in section 3 we point out the 
problem existing in this present statistical method, and in- 
troduce the concept of "generation" of normal pulsars; in 
section 4 we introduce a new statistical method on the spin 
period distribution of normal pulsars; in section 5 we make 
a conclusion. 

2 Usual Statistical Method on the Spin 
Periods of Normal Pulsars at Present 

As is well known pulsar spin periods usually span three 
magnitudes or so from milliseconds to seconds. So peo- 
ple like to do statistics for the log values of the spin peri- 
ods (Manchester, Hobbs, Teoh & Hobbs 2005; Manchester 
2009), or directly under the log coordinate of the spin peri- 
ods (Manchester 2009), because the log function can trans- 
form the different scales of spin periods to the same scale. 
In Fig. Q] we illustrate this method for the normal pulsars 
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Fig. 1 The normal pulsar number distribution ac- 
cording to the log values of spin periods. The data 
are from the ATNF pulsar catalogue for pulsars with 
spin periods larger than 30 ms. See the internet Web: 
http://www.atnf.csiro.au/research/pulsars/psrcat (Manch- 
ester et al. 2005). 

whose spin periods exceed 30 ms (Lyne & Smith 2006) 
in the ATNF pulsar catalogue (Manchester et al. 2005). It 
seems to accord with the following Gaussian distribution 
very well: 

y = yO + ^=e- 2 (—\ (1) 
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where yO = 11.741 ± 3.963, A = 243.163 ± 8.562 
to = 0.678 ± 0.020, xc = -0.225 ± 0.008 are the fit- 
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ting parameters. The coefficient of determination (COD) is 
0.99438. The smaller the relative errors and the closer to 1 
the COD, the better the fitting. Apart from yO, the relative 
errors of the other parameters are all less than 5%. So it 
seems safely to say that the normal pulsars distribute Gaus- 
sianly according to the log values of spin periods. 

3 Concept of "Generations" of Pulsars 

But we should not forget that these normal pulsars were 
not born in the same period. In fact, they are born with 
some birthrate 1^3 per century (Vranesevic et al. 2004; 
Faucher-Giguere & Kaspi 2006) which can be larger five 
times taking into the uncertainties to form the current of 
pulsars in the Galaxy. At the beginning of the current, i.e., 
the newborn pulsars have a distribution of the spin periods 
which nowadays is regarded to be about 0.1-0.5 s for up to 
40% of all pulsars (Vranesevic et al. 2004; Lorimer 2006). 
So maybe we could guess that the distribution of spin peri- 
ods of newborn normal pulsars is pulse form like Gaussian 
distribution with the distribution peak at 0. 1-0.5 s. Certainly, 
the distribution pulse form of spin periods of the newborn 
pulsars could not approach to the infinitely because there 
is a smallest rotation period of pulsars smaller than which 
the pulsar will be broken up and which can approximately 
be estimated as small as 1 .5 ms when the centrifugal force is 
just equal to the gravity at the equator (Lyne & Smith 2006). 

In the initial part of the current of pulsars which can 
be imaged as "initial generation of pulsars", the spin period 
distribution should be like the spin period distribution of 
newborn pulsars. Due to the radiative energy taking out the 
rotation energy of pulsars (Gold 1968, 1969; Pacini 1968), 
the rotation of pulsars will turn slow gradually and this 
pulse-like distribution of spin periods of the initial gener- 
ation will move in the right direction of spin period axis, 
certainly with the possibility that some pulsars could not ra- 
diate the radio pulses owing to the extinction of the polar 
cascade process (Chen & Ruderman 1993; Hibschman, Jo- 
hann & Arons 2001) and the pulse distribution form may 
also change as illustrated in Fig. [2] 

Let us focus on the statistical step sizes in Fig. Q] The 
peak in Fig.[T]is at if) -0 - 225 w 0.6 s approaching to the as- 
sumed 0.1-0.5 s magnitude of the distribution peak of the 
spin periods of the initial generation pulsars. So we can 
guess there are not many generations of pulsars left to the 
peak in Fig.Q] In addition the statistical step size left to the 
peak is even so smaller that it could not contain one pulsar 
generation and in fact to some extent reflect on the distribu- 
tion pulse form inner structure of spin periods of the initial 
several generations of pulsars. Whereas the statistical step 
sizes right to the peak in Fig.[T]are generally very large that 
can contain a few generations in the same statistical step. 
So in fact the statistics right to the peak in Fig. Q]re fleets the 
change rule of pulsar numbers of different generations with 
spin periods. 
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Fig. 2 The sketch map of different pulsar generations at 
different spin periods. 

Now we make clear that the statistical step sizes in Fig.Q] 
are suitable to do different things that reflect different phys- 
ical meanings: one is for the distribution of spin periods 
in the pulse-distribution-form of the initial generations; the 
other is for the distribution of pulsar numbers among dif- 
ferent pulsar generations. The log transformation mixes the 
different study scales of the two different questions, conse- 
quently confuses the two questions. The result is that the 
physical meaning of Gaussian distribution in Fig. [T]is very 
limited. 

4 New Statistical Method on the Spin Periods 
of Normal Pulsars 

Since nowadays the pulsars observed are including various 
generations, we can study the pulsar number distribution 
among different pulsar generations, i.e., we would study the 
pulsar number distribution with the spin periods, not with 
the log values of spin periods. Fig.[3]is the result. It accord 
with the following exponential decay with the very surpris- 
ing COD 0.99992: 

y = Al e~ x/tl + c, (2) 
where Al = 2408.115 ± 15.750, tl = 0.738 ± 0.006 s, 
c = 3.614 ± 1.203 in Fig. [3] The fitting of Fig. [3] is ob- 
viously better than that of Fig. Q] We know for pulsars the 
larger the age, the longer the spin period, and the weaker the 
radiative ability. Fig. [3] tells us that the active pulsar num- 
ber would gradually reduce with the increasing of spin peri- 
ods. Meanwhile it also tells the active pulsar number drops 
with the increasing of time. If the spin period has linear re- 
lationship with the age, then the active pulsar number would 
also drop exponentially with time. But the spin periods do 
not have simple relationship with the time (Lyne & Smith 
2006; Camenzind 2007; Gonthier, Van Guilder & Harding 
2004; Contopoulos & Spitkovsky 2006). In spite of this we 
still can introduce the half-decay spin period for pulsars like 
the half-decay period for radiative elements: the value from 
Fig.[3]is about 0.516 ± 0.004 s. 
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Fig. 3 The pulsar number distribution in different pulsar 
generations. The samples are the same as in Fig.Q] 
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5 Conclusion 

The pulsars and radiative elements belong to very different 
scales of the universe respectively, and have very different 
inner active mechanisms, whereas they have such similar 
decay rules, i.e., their activities drop exponentially though 
with different variables: one with the spin period; the other 
with the time. 

Fig. [3] does not take into count the selective effect of 
observation, the most important selective effect is the lumi- 
nosity, because the radiative capability of pulsars declines 
with age. It is very difficult to count this selective effect. Be- 
sides this, we do not forget that pulsars are the high veloc- 
ity objects (Vranesevic et al. 2004; Arzoumanian, Chernoff 

6 Cordes 2002; Hobbs, Lorimer, Lyne & Kramer 2005), 
which would make the number of pulsars in some region 
of the sky have some change. So Taking into account var- 
ious effects will bring very large uncertainties. In spite of 
all these, we still believe that the law disclosed in Fig. [3] 
is true, partly from the philosophical consideration that the 
different scales in the universe might abide by some similar 
law. 
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